Superstability of the yeast cell cycle dynamics: 
Ensuring causality in the presence of biochemical stochasticity 
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Gene regulatory dynamics is governed by molecular processes and therefore exhibits an inherent 
stochasticity. However, for the survival of an organism it is a strict necessity that this intrinsic 
noise does not prevent robust functioning of the system. It is still an open question how dynamical 
stability is achieved in biological systems despite the omnipresent fluctuations. In this paper we 
investigate the cell-cycle of the budding yeast Saccharomyces cerevisiae as an example of a well- 
studied organism. 

We study a genetic network model of eleven genes that coordinate the cell-cycle dynamics 
using a modeling framework which generalizes the concept of discrete threshold dynamics. By 
allowing for fluctuations in the transcription/translation times, we introduce noise in the model, 
accounting for the effects of biochemical stochasticity. We study the dynamical attractor of the 
cell cycle and find a remarkable robustness against fluctuations of this kind. We identify mecha- 
nisms that ensure reliability in spite of fluctuations: 'Catcher' states and persistence of activity 
levels contribute significantly to the stability of the yeast cell cycle despite the inherent stochasticity. 
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I. INTRODUCTION 

Uncovering the complex mechanisms involved in gene 
regulation remains to be a major challenge. While the 
biochemical processes involved in the expression of a sin- 
gle gene are increasingly well understood, the interplay of 
whole networks of genes poses additional questions and 
it is unclear what level of system-specific detail has to 
be taken into account to describe a gene regulatory net- 
work |2ll l24j . One promising approach to modeling sys- 
tem wide dynamical states of a network is to go to an 
abstract level of description Q which may even include 
discrete deterministic models such as Boolean or thresh- 
old networks 0, |2j| . We here extend a recent model of 
the yeast cell cycle dynamics that is successfully based 
on this approach 01 ■ 

The cell cycle of the budding yeast Saccharomyces 
cerevisiae is a widely studied example of a robust dy- 
namical process [22] • In ^3 > f ne yeast cell cycle was 
modeled in the framework of a discrete threshold net- 
work. From the data in , eleven genes that play a key 
role in the cell cycle process were identified along with 
their known (direct or indirect) interactions. The activity 
of a certain gene is modeled as a two-state system, with 
values 1 (active) or (inactive). Using a threshold model 
of interactions, the biological sequence of activity states 
in the process is exactly reproduced. The authors also 
find considerable dynamical robustness properties that 
can be traced to the properties of the basin of attraction 
of the biological fixed point. 

Remarkably, these results were obtained using a dis- 
crete time model, where each discrete time step t is de- 
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fined by the intervals between activity changes. In the 
model, the activity state of every gene is determined 
solely by the state of its transcription factors at the previ- 
ous time step. It is remarkable, that in this case at least, 
the biochemical stochasticity of gene regulation [Til Il9l | 
can be neglected in the model. 

In particular, as was shown in [Til ], attractors under 
synchronous dynamics can be unstable if stochasticity 
is imposed on the transmission times. In this work, we 
investigate whether the cell-cycle process is stable under 
such perturbations. 

Investigations of dynamical robustness have been dis- 
cussed in a variety of different biolo gica l systems, such as 
segmentation in the fruit fly pi IEL l25j. or two-gene cir- 
cadian oscillators [2(| . Different conceptions of the word 
'robustness' have been used 0. Robustness against mu- 
tations means that a specific process can be performed 
reliably by a system even if some changes to the struc- 
ture of the system are conducted. The yeast cell-cycle is 
remarkably robust in this sense |l7| . Other approaches 
to assessing robustness in biological networks include lo- 
cal stability and bifurcation analyses ||, stability under 
node state perturbation and probabilistic Boolean 

networks [20j. In this work we will concentrate on the 
robustness under stochastically varying processing times 
(for protein concentration buildup and decay) as was con- 
sidered in [l5j . 

Other models of the yeast cell-cycle include molecu- 
lar models of major CDK activities in Start and Finish 
states and of S-phase entrance in Q] . In Q stochastic 
differential equations have been used to fit time-courses 
of protein concentration levels in the yeast cell-cycle net- 
work. 
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II. MODEL DESCRIPTION 

Following ^3] j a network of eleven nodes is used to 
describe the cell cycle process. They are given in table 
0] along with the synchronous sequence of activity states 
recorded in that work. Using a technique introduced in 
we extend that model to include fluctuating trans- 
mission delays and to allow for real numbers for protein 
concentrations levels (0 < c;(i) < 1 for protein i). We 
keep the characteristics of the description of ^3] > that is 
the effect of protein j on the transcription of protein i is 
determined by a discrete activity state ('active' or 'inac- 
tive') of protein j. In our continuous description, we set 
the activity state Si of a protein to 1 if the concentration 
is above a certain threshold (c^(t) > 0.5), otherwise it is 
0. 

The transmission function that determines the tran- 
scription or degradation of protein i is given by 



/»•(*, *«0 = 



o, v, r/,,.s; y (/ 



t d ) > o, 
t d ) < o. 



(1) 



where td is the transmission delay time that comprises the 
time taken by processes such as translation or diffusion 
that cause the concentration buildup of one protein to 
not immediately affect other proteins. The numbers 
determine the effect that protein j has on protein i. An 
activating interaction is described by = 1, inhibition 



by en 



- 1 . If the presence of protein j docs not affect 



expression of protein i, ay = 0. 

If J2j a ijSj{t — td) = 0, the value of ft depends on 
whether the node is modeled as a self-degrader. Self- 
degraders are those nodes that are down-regulated by 
external processes (Cln3, Clnl,2, Swi5, Cdc20/Cdcl4, 
Mcml/SFF). Self-degrader nodes will take a value 
fi{t,td) = whereas the transmission function of non- 
sclf-dcgradcrs is left unchanged, i.e. the last time t when 
fi(t,td) 7^ determines the state at time i. 

We now describe the time evolution of the system of 
genes by the following set of delay differential equations 



dcj(t) 
dt 



fi(t,U) 



Ci(t) 



(2) 



For the simple transmission function given above, this 
equation can be easily solved piecewisc (for every period 
of constant transmission function), leading to charging 
behavior of the concentration levels 



Ci(t > t Q ) = 



1 - (1 - c(to)) exp(-(t 
c(f )exp(-(f - t Q )/r) 



to)M fi>0, 
fi<0. 
(3) 

This has the effect of a low-pass filter, i.e., a signal 
has to sustain for a while to affect the discrete activity 
state. A signal spike, on the other hand, will be filtered 
out. Concentration buildup in our model is depicted in 
figure^ Here, the transcription factor of a protein is as- 
sumed to be present in the time span t = — 1 to (upper 
panel). The production of the protein starts after the 
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FIG. 1: Concentration buildup and decay of a protein given 
a specific input signal and the corresponding activity state 
(t d = 0.8, t = 0.3). 



delay time (here td = 0.8) and the concentration crosses 
the critical level of 0.5 at t = 1 (central panel), switching 
the activity state to "on" (lower panel). In the case of 
very fast build-up and decay (r — ► in cq. J3J) and with 
the delay time set to one (td = 1), we exactly recover 
the synchronous dynamics of |l7j . Thus, our described 
model is a simple generalization of the synchronous case 
to allow for a continuous time description. 

We now ask the following question: Is the original se- 
quence stable under stochastic timing noise (stochasti- 
cally varying signal delay times) or can the noise cause 
the system to assume different states? As the sequence 
from [l^ (reproduced in table [IJ runs into the stationary 
Gi fixed point and an external signal is needed to trig- 
ger the starting state again, we create a repeating cycle 
of states (limit cycle) by explicitly adding the rule that 
Cln3 production is triggered as soon as the final state in 
the synchronous sequence is reached. We will investigate 
whether this limit cycle is inherently stable or whether 
it needs the perfect synchronization of the artificial syn- 
chronous update. 

In this context it is important to note that the sta- 
bility of the complete cell-cycle system also depends on 
the behavior of all other proteins involved. However, the 
stability of the core genes is most important, as they reg- 
ulate the other proteins. Only if the regulators perform 
reliably, the system whole can be robust. 

To compare the time series of our simulations with the 
discrete time steps of the synchronous case, we record 
a time step whenever the system keeps all its activity 
states constant for a time span of at least td/2. With 
every switch of activity states (say, at time to) we check 
whether the transcription of any other protein P is af- 
fected. If so, the concentration level of protein P will 
begin to rise at time to + td + X where \ denotes a uni- 
formly distributed random number between and Xmax 
that perturbs the delay times. 

Our simulation time is not directly related to the ac- 
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TABLE I: The synchronous sequence of states as recorded in Il7jl . We color time steps in which only one switch occurs blue, 
and those with more than one switch red. 



tual time intervals of the biological processes involved. 
However, we are not so much interested in the specifics 
of the time course but rather in the properties of stability 
and for this assessment it is not important how long the 
actual phases take. Our model captures two principles of 
real world gene regulatory networks: Interactions occur 
with a characteristic time delay (denoted by td)', and we 
use continuous concentration levels and implement low 
pass filter behavior due to protein concentration buildup 
with a characteristic time r |lfj |. 



III. RESULTS 

First, we check if the system reproduces the syn- 
chronous sequence under small perturbations of the delay 
time. Thus we stay in the regime where Xmax is sig- 
nificantly smaller than the characteristic protein decay 
or buildup time r. In the main simulation runs we set 
td = 1, t = 0.3 and Xmax = 0.1, but any numbers that 
fulfill Xmax td give the same results. 

We find that the synchronous sequence of states is reli- 
ably reproduced by this stochastic dynamics. Even long 
simulation runs of t < 10 7 cannot push the system out 
of the original attractor. This means that the biological 
sequence is absolutely stable against small perturbations. 

To understand this, we look at the synchronous se- 
quence of states in tablcQ] In steps 2-^3,4-^5,8^9, 
9 -> 10, 11 -> 12, 12 -> 13 (marked blue in the table) 
only a single protein changes its activity state. If all steps 
were of this kind, fluctuations of the event times would 
not be able to destroy the attractor at all. States marked 
in red denote events where multiples switches happen at 
the same time. 

To illustrate this point, let's assume two nodes switch 
their states at times ii and t^ (we call this a 'phase lag'). 
The system thus assumes an intermediate state in the 
time span between t\ and ti. Approximately at time 
ii + td the next switches occur and due to the interme- 
diate state it is possible that proteins switch their states 
which would normally be constant in this step. Because 



of the charging behavior of the concentration levels, these 
'spikes' will be filtered out. The only way to destroy the 
attractor is thus when the phase lags add up in a series 
of steps. This cannot happen in the yeast cycle, however, 
due to the states marked in blue color in the table. When 
only one protein changes its state in a time step, all di- 
vergence of signal times will be reset and the synchrony 
is restored. We therefore call these steps 'catcher states' 
as they remove phase-lags from the system. 

Now that we know that small perturbations cannot 
drive the system out of the synchronous attractor, we 
want to investigate stability under stronger noise. To 
address this question, we have to loosen our definition 
of stability. Up to now, we have requested the system 
to follow the exact sequence of states of the synchronous 
dynamics. It is clear that this strict stability cannot be 
obtained if we increase the noise to be more than half of 
the transmission delay itself, because two nodes switching 
at the same synchronous time step can receive switching 
times that differ by more than id/2. The intermediate 
step taken when only one node has switched obviously 
violates the stability criterion. 

To assess the stability of the system under strong noise, 
we employ a different stability criterion. We let the sys- 
tem run with the sole constraint that the stationary G\ 
state will be assumed regularly for a time span of at least 
td- Any fluctuations occurring inbetween two G\ inci- 
dences will be tolerated, as long as the system finds its 
way to the G\ state of the cell cycle in which growth 
occurs until the cell size signal is triggered. Although 
this might seem too loose a criterion for robust biologi- 
cal functioning, one has to remember that the cell-cycle 
process is also backed up by a system of checkpoints that 
can catch faulty system states. We investigate here the 
inherent stability of the system disregarding these check- 
points but at the same time allowing more variability in 
the sequence. 

Remarkably, with noise of the order of the delay time 
and largely independent of the filter used, the system 
reliably stays in the biological attractor. An example run 
with td = 1, r = 0.3 and \ = 0.9 ran for a time of 10 7 
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FIG. 2: Time course of a run with noise of the order of the 
delay time Xmax = 0.9 td- Black boxes denote active states, 
white means inactive. On a micro-time level the effect of 
fluctuations is visible, but on a larger time scale the dynamics 
is very stable. 



following the biological attractor sequence (in the wider 
sense mentioned above). A typical time span of this run 
is shown in figure [2J This is a surprising result, because 
in general one expects a system to be able to leave its 
attractor sequence under such strong noise if a series of 
multi-switch events (steps 5 to 8) is involved anywhere 
during the sequence. 

Our proposed criterion is not trivially fulfilled: by 
changes in the sequence of switching events or by delaying 
one of several events that occur at the same synchronous 
time step, a new sequence could be triggered. This could 
force the system to jump into one of the other six fixed 
points identified in \u\ without the possibility to return 
to the biological sequence. In figure [21 we show an ex- 
ample of a simulation run with extremely strong noise 
Xmax = 3 td that shows that the system can jump out 
of the attractor. However, it is also apparent that even 
under such strong fluctuations the system runs quite reg- 
ularly until it finally loses its attractor sequence. 

We now quantify the stability of the biological pathway 
under such strong noise. How likely is it for the system 
to lose its biological sequence and to run into a differ- 
ent fixed point? To address this question, we initialize 
the system at the Start state again and check whether it 
completes the cycle. Again, we use the lose criterion de- 
scribed above, which means we only request the system 
to reach the Start state again. In figure 01 we show the 
ratio of erroneous runs of the biological pathway plotted 
against the noise level Xmax- It can be clearly seen that 
for reasonable noise levels the ratio of sequence runs not 
ending in a biological fixed point is very small. In fact, 
even with unrealistically high noise levels of Xmax = 20 
or more (which amounts to arbitrary update times), only 
in a quarter of the runs the system jumps out of the bi- 
ological state sequence. 



FIG. 3: Time course of an example run with strong noise 
Xmax = 3 td- After some repetitions of the biological state se- 
quence the attractor cycle is lost and a fixed point is assumed. 
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FIG. 4: The ratio of runs escaping the biological limit cycle 
plotted against the maximal noise level Xmax- Even for strong 
noise, the fraction of erroneous runs is very small. 



The by far dominating cause for this (very small) in- 
stability is the first step (cf. table|l} where both SBF and 
MBF are activated by Cln3. If the Cln3 concentration 
is degraded before activating the transcription of cither 
SBF or MBF, the system loses the biological sequence. If 
we explicitly force Cln3 activity to sustain long enough to 
make sure that both SBF and MBF are produced, even 
this small instability vanishes and the system assumes 
practically complete stability for all reasonable noise lev- 
els (0.1% erroneous runs at x max = 3t<j)- This super- 
stability is due to the fact that all proteins keep their 
activity states for an extended time. Extremely strong 
noise is therefore needed to delay a single activity switch 
long enough to significantly perturb the system. 

We have tested all results with a wide variety of pa- 
rameters. With a fixed number for the delay time td. 
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only the noise level x max and the characteristic protein 
buildup time r can be adjusted. Our results are com- 
pletely robust against changes of r, even removing the 
filter completely or setting it an order of magnitude larger 
than the delay time does not affect the robustness prop- 
erties described above. 



IV. DISCUSSION 

As we have shown in the previous section, the yeast 
cell-cycle control network is astonishingly stable against 
fluctuations of the protein activation and degradation 
times. The network and the resulting dynamics exhibit 
a number of features that cause this stability: As was 
already discussed in [Tt| , the basin of attraction is very 
large, making it unlikely that an intermediary state be- 
longs to one of the other fixed point basins. A second 
remarkable property is that all node states are sustained 
for at least three (synchronous) steps, making the sys- 
tem less dependent on the specifics of the concentration 
buildup procedure. Third and most important for the 
observed supcrstability under noisy transmission times, 



is the presence of the catcher states which prevent the 
system from gradually running out of synchrony. 

Thus, we have seen that without even taking into ac- 
count the biological checkpoint mechanisms that give ad- 
ditional stability and error-correction features, the sys- 
tem shows a strong inherent robustness against intrin- 
sic fluctuations. In this example of the yeast cell-cycle 
dynamics, potential mechanisms that provide robustness 
under biological noise can be observed. A system with- 
out an external clock (or any other external control) can 
still run reliably if it has intrinsic features that enforce 
robustness: catcher states, persistence of states and an 
attractor landscapes that minimizes the possibilities to 
escape the biological sequence. 

To conclude, we have investigated the stability of the 
cell-cycle network by extending the model of Li et al. to 
allow asynchronous updating of the activity states of the 
genes. We find that the system exhibits robust behav- 
ior under noisy transmission times. Even without taking 
into account the checkpoint mechanisms that give addi- 
tional stability and fallback features, the system shows a 
strong inherent robustness that aids in maintaining reli- 
able functioning. 
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